Temporal Changes in Demography and Genetic Diversity of Oriental Storks at the Stage of Long-Term Captive Propagation and Reintroduction Initiation.

Captive propagation and reintroduction are the major steps in the ex-situ conservation of locally extirpated endangered species in a historical region. In a species restoration project conducted in South Korea, we examined temporal changes in demographics and genetic diversity of oriental storks (Ciconia boyciana). Demographic and genetic data from 1996-2018 were analyzed for 80% of all captive and recently reintroduced individuals. Founder establishment and pair formation induced increases in population size and genetic diversity during the early stage of captive propagation. The degree of genetic diversity was found to become saturated and stable with long-term captive propagation. However, this might be a concern for future genetic diversity of both captive and reintroduced populations simultaneously due to the extraction of captive populations at the early stage of reintroduction. Our findings suggest that periodic evaluation of genetic diversity and selection for releasing individuals, using effective genetic markers, would assist in balancing the genetic diversity of the captive and reintroduced oriental storks at the early stage of reintroduction.

The complete mitochondrial genome of Rana coreana (Anura: Ranidae).

Rana coreana is a brown frog species native to the Korean Peninsula. We characterized the complete mitochondrial genome of the species. The mitochondrial genome sequence of R. coreana is 22,262 bp and comprises 13 protein-coding genes, two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and two control regions (CRs). The CR duplication and gene organization were identical to those observed in Rana kunyuensis and Rana amurensis. A total of 13 protein-coding genes were used to examine the phylogenetic relationships between this species and the genus Rana. R. coreana living on the Korean Peninsula, formed a cluster with R. kunyuensis and R. amurensis, with R. coreana showing the closest phylogenetic affinity for R. kunyuensis.

Development of InDel markers for interspecific hybridization between hill pigeons and feral pigeons based on whole-genome re-sequencing.

Interspecific hybridization occurs among birds, and closely related sister taxa tend to hybridize at a high rate. Genomic hybridization markers are useful for understanding the patterns and processes of hybridization and for conserving endangered species in captivity and the wild. In this study, we developed genomic hybridization markers for the F1 progeny of the sister taxa feral pigeons (Columba livia var. domestica) and endangered hill pigeons (Columba rupestris) (family Columbidae). Using whole-genome re-sequencing data, we performed genome-wide analysis for insertion/deletion (InDel) polymorphisms and validated using primers. We conducted polymerase chain reaction (PCR) and agarose gel electrophoresis to identify species-specific InDels. We produced eight F1 hybrids of hill and feral pigeons, and their samples were tested by re-performing analyses and sequencing using 11 species-specific InDel polymorphisms. Eight InDel markers simultaneously amplified two DNA fragments from all F1 hybrids, and there was no abnormality in the sequencing results. The application of genomic tools to detect hybrids can play a crucial role in the assessment of hybridization frequency in the wild. Moreover, systematic captive propagation efforts with hybrids can help control the population decline of hill pigeons.

Whole-genome sequencing revealed different demographic histories among the Korean endemic hill pigeon (Columba rupestris), rock pigeon (Columba livia var. domestica) and oriental turtle dove (Streptopelia orientalis).

BACKGROUND: The family Columbidae is known as the pigeon family and contains approximately 351 species and 50 genera. Compared to the wealth of biological and genomic information on these Columba livia var. domesteca, information on Columba rupestris and Streptopelia orientalis has been rather limited. The C. rupestris population size is decreasing in Korea. OBJECTIVES: Whole-genome sequencing and identification of population characterization of each species based genome variation on 9 Korean pigeon and dove samples, namely, six hill pigeon (C. rupestris), one rock pigeon (C. livia var. domestica) and two oriental turtle dove (S. orientalis) samples. RESULTS: The whole genome of 9 genotypes were sequenced and mapped to the C. livia reference genome. Sequence alignment showed over 96% identity in C. rupestris and 94% identity in S. orientalis to the reference genome (GenBank assembly accession: GCA_001887795.1). Sequence variations, including single nucleotide polymorphisms (SNPs), insertions and deletions (InDels), and structural variations, revealed that intergenus (Columba vs. Streptopelia) variations were approximately four times higher than intragenus variations (C. livia vs. C. rupestris). Of the two Columba species, C. livia var. domestica is closer to S. orientalis than C. rupestris. Pairwise sequentially Markovian coalescent (PSMC) demographic history analysis revealed that the three species underwent a common population bottleneck between 105 and 120 Kya; since then, the effective population sizes of the rock pigeon and oriental turtle dove have increased. CONCLUSION: The effective population size of the hill pigeon, an Endangered Species of Grade II in Korea, has increased slowly from the second severe bottleneck that occurred approximately 0.5-1.4 × 104 years ago. Our results showed no relationship between copy number variation in the Norrie disease protein (NDP) regulatory regions and plumage color patterns. We report the first comparative analysis of three pigeon genomes.

Phylogenetic relationships between three Korean pit viper Gloydius (Serpentes: Crotalinae) species using mitochondrial DNA genes.

BACKGROUND: Molecular phylogenetic studies of the Asian pit viper genus Gloydius have been widely published in Asia, but Korea population have not been conducted till date. OBJECTIVE: This study aimed to analyze the phylogenetic relationships of three Gloydius species (G. saxatilis, G. brevicaudus, and G. ussuriensis) from Korea with other Gloydius species, based on Cytochrome b and ND4. METHODS: We compared 160 samples representing the three species with those of 17 reference species and their phylogenetic status and genetic diversity were analyzed with concatenated sequences of two mitochondrial DNA. RESULTS: Korean G. brevicaudus and G. saxatilis showed high haplotype diversity and relatively low and moderate nucleotide diversity, respectively. Although G. ussuriensis showed high genetic diversity, it was low in the Baengnyeong Island population. The phylogenetic tree represented two major lineages. One major lineage comprised G. ussuriensis, G. tsushimaensis, G. blomhoffii, and G. brevicaudus. The Chinese G. ussuriensis belonged to the same clade as the Korean G. ussuriensis and was closely related to the Baengnyeong Island population. Moreover, G. tsushimaensis was closely related to G. ussuriensis from southwestern Korean and Jeju Island populations. The other major lineage comprised the remaining 12 species and G. saxatilis. Korean G. saxatilis was closely related to G. saxatilis, G. shedanoensis, and G. intermedius from China. CONCLUSION: The phylogenetic status of the Korean Gloydius species in comparison with the other Gloydius species was identified. We suggesting the conservation management unit for the Baengnyeong Island population, while the current conservation status of Korean G. saxatilis is suggested to be revised to a higher level.

Complete mitochondrial genome of Gloydius saxatilis (Viperidae: Crotalinae) from Korea.

In this study, we sequenced the complete mitochondrial genome of Gloydius saxatilis using Illumina next-generation sequencing. The total length of the mitogenome was 17,223 bp, and contained 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs), two non-coding control regions (CRs), and the origin of light (OL)-strand replication. The genome structure and order of the genes were similar to other Crotalinae species. Phylogenetic analysis based on the 13 concatenated PCGs indicated that G. saxatilis closely related to G. intermedius and, G. shedaoensis.

Yellow sea mediated segregation between North East Asian Dryophytes species.

While comparatively few amphibian species have been described on the North East Asian mainland in the last decades, several species have been the subject of taxonomical debates in relation to the Yellow sea. Here, we sampled Dryophytes sp. treefrogs from the Republic of Korea, the Democratic People's Republic of Korea and the People's Republic of China to clarify the status of this clade around the Yellow sea and determine the impact of sea level change on treefrogs' phylogenetic relationships. Based on genetics, call properties, adult morphology, tadpole morphology and niche modelling, we determined the segregated status species of D. suweonensis and D. immaculatus. We then proceeded to describe a new treefrog species, D. flaviventris sp. nov., from the central lowlands of the Republic of Korea. The new species is geographically segregated from D. suweonensis by the Chilgap mountain range and known to occur only in the area of Buyeo, Nonsan and Iksan in the Republic of Korea. While the Yellow sea is the principal element to the current isolation of the three clades, the paleorivers of the Yellow sea basin are likely to have been the major factor for the divergences within this clade. We recommend conducting rapid conservation assessments as these species are present on very narrow and declining ranges.

Complete mitochondrial genome of two shorebirds (Charadriiformes: Scolopacidae), great knot (Calidris tenuirostris) and bar-tailed godwit (Limosa lapponica).

The mitochondrial genome of Calidris tenuirostris and Limosa lapponica were described using the whole mitochondrial genome obtained from Illumina Next-Generation Sequencing (NGS) technology. Total length of the mitogenome of C. tenuirostris was 16,732bp with slight A+T bias (55.3%). Genome size of L. lapponica was 16,773bp long and A+T biased (56.3%). Both gemones consisting of 2 rRNAs, 13 protein-coding genes, 22 tRNA genes and 1 non-coding regions. This is the first report of complete mitogenomes of these two shorebird species, (C. tenuirostris and of L. lapponica). We observed paraphyletic relationship among the species in the Family Scolopacidae. Also our result showed analogous patterns with the previous studies on the parallel relationships of shorebird species. This study provides basic genetic information for help in understanding phylogenetic relationships . within the Charadriiformes.

A core set of microsatellite loci for yellow-throated marten, Martes flavigula: a case of inferences of family relationships.

BACKGROUND: Microsatellite markers are an ideal molecular marker for population genetic studies such as population structure, pedigree, and kinship. The yellow-throated marten (Martes flavigula) is widely distributed in coniferous and deciduous forests of eastern Asia and plays the role of an indicator and umbrella species in South Korea, given the absence of top predators such as tiger and leopard. OBJECTIVE: The aim of our study was to establish a core set of microsatellite markers that could be used for a population genetics study on M. flavigula. METHODS: We characterized 21 di-motif microsatellites for M. flavigula by Illumina next-generation sequencing. We evaluated them for a population genetics study against five established criteria together with 33 previously developed microsatellites. We calculated relatedness values between individual yellow-throated martens in two groups that were suspected to be siblings using the selected core set of markers to confirm applicability. RESULTS: Twenty-three loci were determined as the core set of microsatellite markers. The probability of identity P(ID) and probability of identity between siblings P(ID)sib of the core set was estimated as 2-15 and 2.2-7, respectively. Relatedness values between individuals in the two groups of M. flavigula revealed that one of the pairs was sisters, confirming that the core set can be applied to kinship studies. CONCLUSION: The developed microsatellite core set in this study is expected to contribute to studies on molecular ecology and population structure of M. flavigula.

The complete mitochondrial genome information of Phoxinus phoxinus (Cypriniformes: Cyprinidae) on the Korean Peninsula and the phylogenetic implication.

Phoxinus phoxinus is a small Leuciscinae species predominantly found in cool and well-oxygenated streams throughout a wide area encompassing Europe, Siberia and East Asia. It is believed that the populations in Korea hold important clues to how the species has been distributed south along the Eurasian continent to the Korean Peninsula. We characterized the complete mitochondrial genomes of two individual fin-clip samples collected from the two Korean river systems. The whole sequences were 17,665 and 18,220 bp, respectively, and included 13 protein-coding genes, 2 ribosomal RNA genes and 22 transfer RNA genes. The genome size difference was due to the considerably different sizes of the control region. The overall genome structures were identical to those observed in other Leuciscinae species.

Microsatellite DNA analysis reveals lower than expected genetic diversity in the threatened leopard cat (Prionailurus bengalensis) in South Korea.

To optimize conservation efforts, it is necessary to determine the risk of extinction by collecting reliable population information for a given species. We developed eight novel, polymorphic microsatellite markers and used these markers in conjunction with twelve existing markers to measure genetic diversity of South Korean populations of leopard cat (Prionailurus bengalensis), a species for which population size and habitat area data are unknown in the country, to assess its conservation status. The average number of alleles and the observed heterozygosity of the species were 3.8 and 0.41, respectively, and microsatellite diversity was lower than the average genetic diversity of 57 populations of 12 other felid species, and lower than that of other mammal populations occurring in South Korea, including the raccoon dog (Nyctereutes procyonoides), water deer (Hydropotes inermis), and endangered long-tailed goral (Naemorhedus caudatus). Furthermore, analysis of genetic structure in the national leopard cat population showed no clear genetic differentiation, suggesting that it is not necessary to divide the South Korean leopard cat population into multiple management units for the purposes of conservation. These results indicate that the genetic diversity of the leopard cat in South Korea is unexpectedly low, and that the risk of local extinction is, as a result, substantial. Thus, it is necessary to begin appropriate conservation efforts at a national level to conserve the leopard cat population in South Korea.

Complete mitochondrial genome sequence of Macromia daimoji Okumura, 1949 (Odonata: Macromiidae).

The dragonfly Macromia daimoji Okumura, 1949 (Odonata: Macromiidae) has been listed as an Endangered insect in South Korea. We sequenced the complete 15,198 bp mitochondrial genome (mitogenome) of this organism, which is the first mitogenome sequence reported from the family Macromiidae. The genome includes a typical set of genes [13 protein-coding genes (PCGs), 2 rRNA genes, and 22 tRNA genes) and one non-coding region with an arrangement identical to that observed in most insect genomes. Phylogenetic analyses using concatenated sequences of the 13 PCGs and 2 rRNA genes using the Bayesian inference (BI) method placed Macromiidae, represented by M. daimoji, as a sister group to Libellulidae with the highest nodal support [Bayesian posterior probabilities (BPP) = 1]. Unlike conventional phylogenetic analysis, the suborders Anisozygoptera and Zygoptera formed a strong sister group (BPP =1), justifying the use of different molecular markers for phylogenetic analysis.

Complete mitochondrial genome of the Northern Long-eared Owl (Asio otus Linnaeus, 1758) determined using next-generation sequencing.

In this study, the mitogenome of Asio otus, the Northern Long-eared Owl, was analysed using Illumina next-generation sequencing. The mitogenome was found to be a circular molecule, 17,735 bp long with a slight AT bias (53.0%). The gene arrangement pattern was the same as that of a typical vertebrate, containing 37 genes (13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and a non-coding control region). In the putative control region (1984 bp), there were two types of tandem repeats at the end of the region. A phylogenetic tree was constructed using the 13 PCG sequences discovered in this study and those of that have been previously published of other Strigidae species and revealed a close relationship between A. otus and A. flammeus. The newly generated mitogenome from this study enriches the genomic resources available for future evolutionary studies and promotes conservation genetics of this species.

Complete mitochondrial genome and phylogenetic analysis of Callorhinus ursinus: an endangered species from South Korea.

We present the complete mitochondrial genome and a phylogenetic analysis of Callorhinus ursinus, the northern fur seal, determined using Illumina next-generation sequencing (NGS) technology. The total length of the mitogenome was 17,154 bp, which consisted of 13 protein-coding genes, two ribosomal RNA genes, 22 tRNA genes, and one control region. The base composition of the entire mitogenome was 33.5% (A), 26.3% (C), 13.9% (G), and 26.3% (T) with an A + T bias of 59.8%. The control region contained two types of tandem repeats. A neighbour-joining (NJ) tree was constructed and comprised two clades with C. ursinus forming a monophyletic group. Data produced in this study will aid exploration of the genetic diversity of endangered C. ursinus and contribute to molecular identification of this species.

The mitochondrial genome of Milvus migrans (Aves, Accipitriformes, Accipitridae), an endangered species from South Korea.

The complete mitochondrial genome sequence of black kite, Milvus migrans, one of the most common diurnal raptor, was characterized using next generation sequencing. The whole genome size was 18,016 bp and consisted of 13 protein-coding genes, 22 tRNAs, 2 rRNAs, a putative control region (CR), and a second control region (pseudo-CR). A frameshift mutation was found in the ND3 gene. Phylogenetic analysis illustrated monophyly of the subfamily Melieraxinae with high statistical support. The genetic resource obtained here will help to resolve taxonomic issues related to subspecies of M. migrans and will act as a starting point for conservation genetics.

Myotis rufoniger genome sequence and analyses: M. rufoniger's genomic feature and the decreasing effective population size of Myotis bats.

Myotis rufoniger is a vesper bat in the genus Myotis. Here we report the whole genome sequence and analyses of the M. rufoniger. We generated 124 Gb of short-read DNA sequences with an estimated genome size of 1.88 Gb at a sequencing depth of 66× fold. The sequences were aligned to M. brandtii bat reference genome at a mapping rate of 96.50% covering 95.71% coding sequence region at 10× coverage. The divergence time of Myotis bat family is estimated to be 11.5 million years, and the divergence time between M. rufoniger and its closest species M. davidii is estimated to be 10.4 million years. We found 1,239 function-altering M. rufoniger specific amino acid sequences from 929 genes compared to other Myotis bat and mammalian genomes. The functional enrichment test of the 929 genes detected amino acid changes in melanin associated DCT, SLC45A2, TYRP1, and OCA2 genes possibly responsible for the M. rufoniger's red fur color and a general coloration in Myotis. N6AMT1 gene, associated with arsenic resistance, showed a high degree of function alteration in M. rufoniger. We further confirmed that the M. rufoniger also has bat-specific sequences within FSHB, GHR, IGF1R, TP53, MDM2, SLC45A2, RGS7BP, RHO, OPN1SW, and CNGB3 genes that have already been published to be related to bat's reproduction, lifespan, flight, low vision, and echolocation. Additionally, our demographic history analysis found that the effective population size of Myotis clade has been consistently decreasing since ~30k years ago. M. rufoniger's effective population size was the lowest in Myotis bats, confirming its relatively low genetic diversity.

The mitochondrial genome of the long-billed plover, Charadrius placidus (Charadriiformes: Charadriidae).

The present study provides the first full mitochondrial genome sequence of a long-billed plover from South Korea. This mitogenome is 16,895 bp in length and consists of 13 protein-coding genes, 2 ribosomal-RNA genes, 22 transfer-RNA genes, and a non-coding control region. Overall, base composition was: A, 31.4%; C, 31.1%; G, 13.8%; and T, 23.8%. A frameshift mutation in the ND3 gene was identified and a tandem repeat (AACA) was observed in the D-loop region. The phylogenetic analysis based on concatenated coding genes indicated Charadriidae's monophyly. These results contribute to further understanding of long-billed plover phylogenetic relationships and species identification.

Complete mitochondrial genome of Haematopus ostralegus (Charadriiformes: Haematopodidae).

The Eurasian oystercatcher (Haematopus ostralegus), Near Threatened in the IUCN red list, was designated classified endangered species II by the Ministry of Environment of Korea and a Natural Monument (No. 326) by the Cultural Heritage Administration of Korea. In this study, the complete mitochondrial genome (16,798 bp) of H. ostralegus was determined for the first time, including 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and 1 non-coding region. The overall base composition was A (31.4%), C (31.3%), G (13.8%), and T (23.5%), so the percentage of A and T (54.9%) was slightly higher than that of G and C. A phylogenetic analysis using concatenated mitogenomes revealed that the family Haematopodidae has a closer relationship with the family Recurvirostridae and H. ostralegus and H. ater clustered together. The results are expected to provide useful resources for species identification and further phylogenetic studies of genus Haematopus.

Sequencing and analysis of the complete mitochondrial genome of Hyla suweonensis (Anura: Hylidae).

The present study reports the complete mitochondrial genome of the Suwon tree frog Hyla suweonensis from South Korea. This endangered species is endemic to Korea. The circular mitogenome of H. suweonensis includes 16,895 bp length and contains 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a non-coding region, which is the typical gene arrangement found in the available Hylinae mitogenomes. Phylogenetic analysis of six Hyla spp. mitochondrial genomes revealed that H. suweonensis is closer to H. tsinlingensis.